ORBITAL ANGULAR MOMENTUM PARTON 
DISTRIBUTIONS IN QUARK MODELS 



SERGIO SCOPETTA AND VICENTE VENTO 

Departament de Fisica Tedrica, Universitat de Valencia, 46100 Burjassot 

(Valencia), Spain 
E-mail: sergio.scopetta@uv.es, vicente.vento@uv.es 

At the low energy, hadronic, scale we calculate Orbital Angular Momentum (OAM) 
twist-two parton distributions for the relativistic MIT bag model and for non- 
relativistic quark models. We reach the scale of the data by leading order evolution 
in perturbative QCD. We confirm that the contribution of quarks and gluons OAM 
to the nucleon spin grows with Q 2 , and it can be relevant at the experimental scale, 
even if it is negligible at the hadronic scale, irrespective of the model used. The 
sign and shape of the quark OAM distribution at high Q 2 may depend strongly on 
the relative size of the OAM and spin distributions at the hadronic scale. Sizeable 
quark OAM distributions at the hadronic scale, as proposed by several authors, 
can produce the dominant contribution to the nucleon spin at high Q 2 . 

1 Introduction 

Understanding how the partons carry the angular momentum in the nucleon 
has become a main effort of present day physics. The quark spin contribution 
AE is well defined in QCD u, thus measurable; with respect to the gluon 
spin contribution Ag, the experimental □ and theoretical □ situation has been 
discussed at the Conference. Our knowledge of the quark, L q and gluon, L g , 
Orbital Angular Momentum (OAM) is less satisfactory. 

It is well known, that the most natural definition of OAM for quarks and 
gluons cannot be separated in a gauge invariant way from the correspond- 
ing spin terms □. However, recently, new definitions of angular momentum 
have been implemented to accommodate gauge invariance EM, and from them 
adequate twist two OAM distributions have been constructed. 

In this respect, three different possibilities have been investigated. One 
proceeds by choosing a particular gauge □'[}□. The OAM operator leads to the 
naive definition of orbital angular momentum up to effects not controlled by 
the gauge fixing. A second, maintains gauge invariance, by loosing covariance, 
defining the distributions only in the class of reference frame where the nucleon 
has a definite polarization In this case the resulting distributions can. lip 
related to the forward limit of off- forward quantities, and are measurable Ehj EI. 
The last proceeds by defining OAM operators such that the distributions are 
gauge invariant □ . Furthermore, in the light-cone gauge, they reduce to 
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the natural definitions li. At present, however, no physical process has been 
proposed to access them. 

hition equations for the OAM distributions have been derived 
and they have been numerically solved by using as input data- 
inspired OAM distributions lia. 

During the last few years we .have developed a scheme to study distribu- 
tions based on model calculations lL-3. This procedure has been worked out thus 
far to leading order in the twist expansion and therefore we should compare 
with the data only at high Q 2 , where the contribution from the non lead- 
ing twists vanishes. In this talk we present the OAM distributions obtained 
by evolving those, properly calculated within relativistic and non relativistic 
models, to the experimental scale Ej. 



2 The theoretical framework 

It has been suggested iiJ-the past that the OAM contribution to the nucleon 
spin can be large at lowE30, as well as large EJ, energy scales. Since the quark 
spin, AE(Q 2 ), and the gluon spin, Ag(Q 2 ), are observables, one may use the 
Spin Sum Rule □ for the nucleon, 



-AE(Q 2 ) + Ag(Q 2 ) + L q (Q 2 ) + L g (Q 2 ) = - (1) 

to determine the global OAM contribution to the spin. The problem of sep- 
arating the OAM in the quark and gluon fractions in a gauge invariant way, 
already addressed in □, is a cumbersome one. 

We will adopt here the gauge independent, twist-2, new definition for 
the quark OAM distribution, given in ref. This definition mixes the po- 
larized and unpolarized singlet quark distributions with the forward limit of 
off- forward parton distributions (OFFPD). Since the OAM matrix element is 
not a Lorentz scalar, there is an ambiguity in this definition for any relativistic 
quantum theory To avoid it, one must take for QCD a system of coordinates 
where the nucleon has a definite helicity, thus loosing covariance. 

We proceed to use this last theoretical development, to perform a phe- 
nomenological model analysis of the OAM distributions, using relativistic as 
well as non-relativistic quark models. 

Let us first discuss the definition of the OAM for a relativistic model. 
According to Ref. if we assume that the nucleon is moving in the z direction 
and is polarized with helicity +1/2 , the quark OAM distribution is given by, 
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L q (x, Q 2 ) = - [x(H(x, Q 2 ) + E q (x, Q 2 )) - A£(z, Q 2 )] , (2) 

where Ti{x, Q 2 ) (AE(i, Q 2 )) is the usual unpolarized (polarized) singlet quark 
distribution and E q (x, Q 2 ) is the f onward limit of the helicity-flip, chiral odd, 
twist-two OFFPD E(x, A 2 , A ■ n) U. The latter quantity is defined through 
the twist-2 part in the the following twist expansion 



/ 



^(P'|V(-An/2) 7^ ^(An/2)|P) = H(x, A 2 , A • n)U(P')^U(P) 

+ E(x,A 2 ,A-n)U(P') l -^^U(P) + ..., (3) 

where P(P') is the 4-momentum of the initial (final) nucleon in a virtual 
Compton scattering process, A M = P' — P M , n = (1, 0, 0, — 1)/2A, A is fixed by 
the choice of the reference frame, and H(x, A 2 , A-n) is the helicity-conserving, 
chiral even, twist-2 OFFPD whose forward limit is the usual forward unpo- 
larized parton distribution. In the definitions above, the forward limit corre- 
sponds to A 2 -> 0, A • n -> 0. 

In the non-relativistic case, the nucleon wave function is given in general 
by an expansion in terms of the eigenstates of some approximate hamiltonian. 
Let N label the quantum numbers of the eigenstates, i.e., principal quantum 
number, orbital angular momentum, spin, then 

*(pi,P2,P3) = 22a N if)N(pi,p3,p3) (4) 

N 

The OAM quark parton distribution, generalizing our approach developed for 
unpolarized and polarized distributions La, is determined by Ej 

L q (x,Q 2 ) = 2irMj2\aN\ 2 / dppL z qtN (p) , (5) 

N J\p-{x)\ 

where 



L q,N(p) = (4>N(pi,P2,P3)\^L*6{p-pi)\-ip N (pi,P2,P3)), (6) 

i=l 

and the lower integration limit is fixed by energy conservation Q. 
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Having set down the framework which defines the OAM distributions at 
the hadronic scale, we proceed to calculate them explicitly in two models and 
to study their evolution. 

3 OAM parton distributions at high Q 2 

We proceed to study the OAM distribution in two different scenarios for pro- 
ton structure: i) a non-relativistic scheme based on the Isgur-Karl model t3; 
ii) a relativistic approach, as described by the MIT bag model Ej. 

For the non-relativistic scenario we consider initially the Isgur-Karl model 
with a proton wave function given by a harmonic oscillator potential including 
contributions up to the 2Two shell E3. In this case the wave function, Eq. (^), 
is given by 

\N) =a s \ 2 S 1/2 )s + a s ,\ 2 S' 1/2 )s + a M \ 2 S 1/2 ) M + a v \ 4 D 1/2 ) M , (7) 

where we have used the spectroscopic notation \ 2S+1 Xj)t, with t = A,M,S 
being the symmetry type. The coefficients were determined by spectroscopic 
properties to be: as — 0.931, as> = —0.274, a» = —0.233, a-p — —0.067. 

Calculating Eq. (^) usiag the wave function Eq. (Q), one gets for the 
parton distribution, Eq. (JfpLj, 

M^) = K| 2 ^ 5 ) (w + 30^ + 45j e ^' (8) 

where a 2 is a parameter of the model . Note that only the small \ 4 D 1 / 2 )m 
wave component gives a contribution to the distribution. 

In the relativistic scenario, we evaluate the twist-two distribution Eq. (^J) 
in the MIT—bag model. The term A£(x, /Zq) in Eq. (Q) is discussed for the 
bag in ref. El To calculate Eq.(||) we need the quantity S(a;, fx%) + E q (x, 
which is given by the forward limit of the OFFPD distribution H(x, A 2 , A • 
n) + E(x, A 2 , A • n). The result for the latter in the bag, calculated by using 
the general definition (Eq. (j|)), is to be found in Eq.(29) of ref. cil, from 
which we have obtained thp_forward limit. The explicit expression for such a 
limit can be found in Ref. Ii3. 

In both cases, Isgur-Karl (hK) and MIT, we use the corresponding support 
correction as defined in tfl and E3, respectively. 

Once the distributions L q (x, /i 2 ,) have been obtained at the low (hadronic) 
scale of the model, w£_perfor_m a LO QCD evolution according to the equations 
displayed in refs. BHlld'tj'EJ. These contain a complicate mixing between 
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Figure 1. Proton OAM distributions in the IK model for the quarks, (a), and for the 
gluons, (b). The full curves show the initial distributions at the hadronic scale of /ig = 0.08 
GeV 2 , where a negligible fraction of the nucleon momentum is carried by the gluons; the 
dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the long-dashed 
curves give the LO-evolved distributions at Q 2 = 1000 GeV 2 . 



Lq(g) ( x > A*o)> A£(x,/i§) and /S.g(x,^), a feature which will be very relevant 
in the analysis of the data. 

The results of our analysis are shown in Figs. 1 through 5. In Fig. 1 we 
show the IK result for quarks (gluons) in (a) ((b)); The full curve corresponds 
to the initial distribution, which is missing in (b), since we start from a very 
low hadronic scale where no valence gluons exist. At LO the hadronic scale 
corresponds to /z§ ~ 0.08 GeV 2 . The dashed curves correspond to the result 
of the evolution from the hadronic OAM distributions, to Q 2 = 10 GeV 2 
(short-dashed) and to Q 2 = 1000 GeV 2 (long-dashed). 

We can summarize the results of the calculation as follows: i) the evolved 
distributions are negative; ii) the magnitude of the distributions increases with 
Q 2 at low x; hi) the magnitude of L q (x, Q 2 ) is small but increases with respect 
to the tiny starting distribution (Eq. (^J)); iv) though gluons are assumed 
to be negligible at the hadronic scale, L g (x,Q 2 ) becomes much larger than 
L q (x,Q 2 ) at high Q 2 . We show in Fig. 2 the same analysis for the MIT bag 
model. Here the initial L q is much larger, but the result of the evolution is 
qualitatively basically the same. 

Thus our hrst conclusion is that there is little model dependence for dif- 
ferent initial OAM distributions. In this case their structure does not seem 
to influence very much the evolution. It is clear that the other inputs of the 
equations, A£(x,/!q) and Ag(x,/x§), are the dominating features. 
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Figure 2. The same as in Fig.l, but for the MIT model. It should be noticed that the 
distribution at the hadronic scale of fj,Q = 0.08 GeV 2 is shown as it comes out from the 
model calculation, before the support correction is implemented. For this reason it does 
not go to zero at x = 1. 



For these results the initial scale has been very low and therefore no initial 
gluon distribution was required. In Fig. 3 we show the result of the evolution 
for the IK model starting from a higher scale, /j,q ~ 0.23 GeV 2 . At this 
hadronic scale, for LO, about 40 % of the proton momentum must be carried 
by the gluons. The polarized gluon contribution is built starting from the 
valence quark distributions as done in Ell and suggested in Ell, and we define 
L g (x, /Xq) from L q (x,^) using the same prescription. 

Again the same features as in the first analysis are found. Recapitulating, 
our results show that the input distributions L q (x, /Xq) and L g (x, /Xq) do not 
seem to determine the behavior of their evolved ones, which turns out to be 
governed by the singlet polarized distributions Ag(x, (Aq) and AE(x,/Xg), due 
to their mixing in the evolution equations. To check to what extent such a 
statement, stressed also in 113, is valid, we analyze in Fig. 4 the results of 
.a_modification of the IK model, the so called D model, already studied in 
cB. In this variant model, the D wave probability is set large to reproduce 
the axial coupling constant of the nucleoli This condition requires the 
following choice of the parameters in Eq. (^): as = 0.894, as> — 0, clm = 
0, aj) — —0.447, i.e., the probability to find a nucleon in the D wave is 
about 20 %. Moreover we have introduced also in the D Model scenario 
polarized valence gluons, as we did for the IK scenario of Fig. 3. From the 
figure it is clear that, while the result for the gluons does not differ in a 
relevant way from the ones found before with the various models, the result for 



ts99: submitted to World Scientific on February 1, 2008 6 



Figure 3. Proton OAM distributions in a modified IK model (see text) for the quarks, (a), 
and the gluons, (b). The full curves show the initial distributions at the hadronic scale of 
[1q = 0.23 GeV 2 , where around 40 % of the nucleon momentum is carried by the gluons; the 
dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the long-dashed 
curves give the LO— evolved distributions at Q 2 = 1000 GeV 2 . The initial gluon distribution 
is so small, in relation with the final one, that it does not show up in (b). 



L q (x,Q 2 ) does, this distribution becoming rather large and positive for large 
x. It is important to stress that in order for this to occur two mechanisms 
were needed, a large initial OAM distribution (as in the MIT bag Model), 
and a higher hadronic scale (as provided by the valence gluons). As shown 
previously, the independent action of the two mechanisms does not lead to 
this behavior. Thus, /S.g(x,iJ^) and A£(:e,/Xq) are governing the evolution as 
long as they are much larger than L q (x,/j,Q) at the initial scale. When they 
have similar size the above statement is not true any more. Note that the IK 
interaction together with the choice of parameters of the D model does not 
describe the hadron spectrum. Nonetheless, other models of interaction (for 
example E130) predict at the low scale a large OAM and fit the spectrum. 
We conclude therefore that a precise knowledge of the OAM distributions will 
serve to distinguish among the models. In Fig. 5 we show the evolution of the 
various contributions to the spin sum rule, Eq. (Q). Fig. 5 (a) corresponds to 
the modified IK scenario used in Fig. 3, where at the scale of the model the 
OAM carried by quarks and gluons is very small and the rest of the proton spin 
is almost equally shared between quarks and gluons spins (L 9 (/Xg) + L g (p^) ~ 
0.01, AE(jug) ~ 0.48 and Ag(nl) ~ 0.25), whereas Fig. 5 (b) corresponds 
to the extreme scenario used already in Fig 4, the so called D model, with 
a large initial OAM (L q (jj%) = 0.145, L g (pl) = 0.055, AS(/ig) ^ 0.4 and 
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Figure 4. Proton OAM distributions for the "D-Model" (see text), for the quarks, (a), 
and the gluons, (b). The full curves show the initial distributions at the hadronic scale of 
fiQ = 0.23 GeV 2 , where around 40 % of the nucleon momentum is carried by the gluons; the 
dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the long-dashed 
curves give the LO-evolved distributions at Q 2 = 1000 GeV 2 . 



A<7(//q) ~ 0.1). As predicted by total angular momentum conservation U 
and already obtained in as a model-independent feature of the evolution 
equation , it is seen that at large Q 2 the huge negative contribution L g (Q 2 ) 
basically cancels out with the positive Ag(Q 2 ). Anyway, the role of the quark 
OAM is found to be very important in the second scenario (cf. Fig. 5(b)), 
being at large Q 2 the largest contribution to the saturation of the spin sum 
rule. Again, we see that quark OAM, due to evolution, can be important at 
large Q 2 if it is not negligible at the scale of the model, whereas the gluons 
OAM, though it is large, is basically cancelled by the gluons spin. 

4 Conclusions 

We have studied the OAM distributions as defined newly in order to take into 
account gauge invariance. _ 

We have seen that evolution, as in previous calculationsliJ, plays a major 
role in the outcome of the predictions. The fact that the gluon and quark 
spin singlet distributions mix in these equations with the quark OAM distri- 
butions, implies that for large Q 2 large contributions from the OAM are to 
be expected, even if they are not present at the hadronic scale. Thus, two 
scenarios arise in a natural way. One, the more conventional one, as described 
by the more traditional models, is defined by quark OAM distributions at a 
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Figure 5. The contributions to the proton spin sum rule, Eq. (gj), according to: (a) the 
modified IK scenario of Fig. 3; (b) the "D model" scenario of Fig. 4. The dashed curve 
shows iAS(Q 2 ), the long-dashed one Ag(Q 2 ), the dot-dashed curve is L q (Q 2 ), the dot- 
long-dashed curve gives L g (Q 2 ) and the full curve represents the sum of the previous four 
terms, giving the spin sum rule (J = i). 



small liadronic scale. In this scenario the evolved distributions are large, neg- 
ative and almost model independent and angular momentum DIS physics is 
dominated by the quark and gluon spin singlet distributions, not by OAM 
distributions at the hadronic scale. A second scenario is defined by quark 
and gluon OAM distributions at a higher hadronic scale. In the latter, soft 
evolution scenario, the initial distributions are important and therefore DIS 
physics may be able to discriminate between models. If the OAM distribu- 
tions are large the outcome of the evolution is strongly dependent on the 
initial distributions and completely different from that of the first scenario. 
Finally the gluon OAM distributions become huge through evolution, even 
they are not present at the hadronic scale. However, as it is well known 
, the gluon OAM and gluon spin contributions cancel to a great extent in 
the nucleon spin, but not so in other moments. Our past experience suggests 
that LO results provide a reasonable qualitative approximation and we do not 
expect that NLO corrections can spoil their general features. 

Many phenomenologicaLimplications have arisen of our study. A careful 
analysis of gauge invariance HEm has permitted us to obtain many observables, 
which may not only lead to a better understanding of the proton spin, but 
to describe the proper behavior of QCD at low energies, i.e., in the confining 
region. These observations are instrumental in defining the picture of the 
proton that should be used for describing low energy properties. 
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